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A NEW DIRECTION IN HYDRODYNAMIC STABILITY: 

BEYOND EIGENVALUES 


Lloyd N. Trefethen 1 , Anne E. Trefethen 2 , Satish C. Reddy* , Tobin A. Driscoll 4 


ABSTRACT 

Fluid flows that are smooth at low speeds become unstable and then turbulent at higher 
speeds. This phenomenon has traditionally been investigated by linearizing the equations 
of flow and looking for unstable eigenvalues of the linearized problem, but the results agree 
poorly in many cases with experiments. Nevertheless, it has become clear in recent years 
that linear effects play a central role in hydrodynamic instability. A reconciliation of these 
findings with the traditional analysis can be obtained by considering the “pseudospectra” of 
the linearized problem, which reveal that small perturbations to the smooth flow in the form 
of streamwise vortices may be amplified by factors on the order of 10 5 by a linear mechanism, 
even though all the eigenmodes are stable. The same principles apply also to other problems 
in the mathematical sciences that involve non-orthogonal eigenfunctions. 
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Introduction 

A change of paradigm is taking place in the field of hydrodynamic stabil- 
ity. This is the field devoted to the study of how laminar fluid flows become 
unstable, the precursor to turbulence (1, 2). It is well known that turbulence 
is an unsolved problem, but not so well known that despite the efforts of gen- 
erations of applied mathematicians beginning with Stokes, Kelvin, Rayleigh, 
and Reynolds a century ago, many of the presumably simpler phenomena of 
hydrodynamic stability also remain ill-understood (3). 

The traditional paradigm is eigenvalue analysis: 

I. Linearize about the laminar solution; 

II. Look for unstable eigenvalues of the linearized problem. 

An “unstable eigenvalue” means an eigenvalue in the complex upper half-plane, 
corresponding to an eigenmode of the linearized problem that grows exponen- 
tially as a function of time t. It is natural to expect that a flow will behave 
unstably if and only if there exists such a growing eigenmode, and over the 
years a great deal of knowledge has been obtained about which flows possess 
such modes, a distinction that depends on the geometry, the Reynolds number, 
and sometimes other parameters. 

For some flows, the predictions of eigenvalue analysis match laboratory ex- 
periments. Examples are Rayleigh- Ben ard convection (a stationary fluid heated 
from below) and Taylor-Couette flow (between a stationary outer and a rotating 
inner cylinder). For other flows, the predictions of eigenvalue analysis fail to 
match most experiments. In this article we consider the two most extensively 
studied examples of this kind: (plane) Couette flow, the flow with a linear ve- 
locity profile between two infinite flat plates moving parallel to one another, 
and (plane) Poiseuille flow, the flow with a parabolic velocity profile between 
two stationary plates (Fig. 1). Other examples where eigenvalue analysis fails 
include pipe Poiseuille flow (in a cylindrical pipe) and, to a lesser degree, Blasius 
boundary layer flow (near a flat wall). 

For Poiseuille flow, eigenvalue analysis predicts a critical Reynolds number 
(4) R = 5772.22 at which instability should first occur (5), but in the laboratory, 
transition to turbulence is observed at Reynolds numbers as low as 1000 (6). 
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Fig. 1 . Velocity profiles for the two laminar flows considered in this article 
(independent of x and z). The geometry is an infinite 3D slab of viscous 
incompressible fluid bounded by parallel walls. The laminar solutions sat- 
isfy the Navier-Stokes equations for all Reynolds numbers R, but for higher 
R the flows are unstable and rapidly become turbulent. 


For Couette flow, eigenvalue analysis predicts stability for all R, but transition 
is observed for Reynolds numbers as low as R ss 350 (7). These anomalies of 
“subcritical transition to turbulence” have been recognized for many years, and 
the explanation has usually been attributed to step I above. If linearization has 
failed, the reasoning goes, one must look more closely at the nonlinear terms, 
or perhaps linearize about a solution other than the laminar one (the theory of 
“secondary instability” (8, 9)). 

Recently it has emerged, however, that the failure of eigenvalue analysis 
may more justly be attributable to step II. It is a fact of linear algebra that 
even if all of the eigenvalues of a linear system are distinct and lie well inside 
the lower half-plane, inputs to that system may be amplified by arbitrarily large 
factors if the eigenfunctions are not orthogonal. A matrix or operator whose 
eigenfunctions are orthogonal is said to be normal (10, 11). The linear operators 
that arise in the Benard and Taylor-Couette problems are normal or close to 
normal. By contrast Reddy, Schmid and Henningson discovered in 1990 that the 
operators that arise in Poiseuille and Couette flow are in a sense exponentially 
far from normal (12). At about the same time the startling discovery was made 
by Gustavsson (13) and Henningson (14) and independently by Butler and 
Farrell (15) that small perturbations to these flows may be amplified by factors 
of many thousands, even when all the eigenvalues are in the lower half-plane. 
The paper by Butler and Farrell is an important and elegant one. It discusses 
many details omitted here, and together with a related more recent paper by 
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Fig. 2. Maximum possible amplification of 3D perturbations in linearized 
Poiseuiile flow as a function of Reynolds number R and x-z wave number 
magnitude k = y/a 2 + (3 2 . In the black region, with leftmost point R = 5772, 
unstable eigenmodes exist and unbounded amplification is possible. Ihe 
contours outside that region correspond to finite amplification factors of 
10 3 , 10 4 (dashed), and 1 x 10 5 ,2 x 10 5 ,...,1.3 x 10 6 . For example, ampli- 
fication by a factor of 1000 is possible for all R > 549. In the laboratory, 
transition to turbulence is observed at R ~ 1000. The analogous picture 
for Couette flow looks qualitatively similar except that there is no black 
region. 


Reddy and Henningson (16), it forms the foundation of the present work (17). 

Fig. 2 illustrates these factors of thousands. The black region has appeared 
in many papers and books and corresponds to parameters for which unstable 
eigenmodes exist (18). The contours defining the shaded regions are new and 
quantify the non-modal amplifications that, may occur in these flows. (A precise 
explanation of this plot is given in the next section.) The possibility of ampli- 
fication of perturbations of viscous flows by non-modal linear mechanisms has 
been recognized for a century (19 23), but until the recent developments, it was 
not known that the magnitudes involved were huge. 
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An essential feature of the amplification process of Fig. 2 is that it acts 
upon 3D perturbations of the laminar flow field. In much of the theoretical 
literature of hydrodynamic stability, attention has been restricted to 2D ( x-y ) 
perturbations, and in particular, the well-known Orr-Sommerfeld equation is 
an eigenvalue equation for 2D perturbations. A principal justification for this 
restriction has been Squire’s theorem of 1933, which asserts that if a flow has an 
unstable 3D eigenmode for some R, then it has an unstable 2D eigenmode for 
some lower value of R (24). The new results indicate that this emphasis on 2D 
perturbations has been misplaced. When only 2D perturbations are considered, 
amplification is still possible, but far weaker. 

This emphasis on 3D, linear but non-modal effects is the emerging new 
paradigm in hydrodynamic stability. 


Streamwise vortices and streaks 

The flow features associated with this amplification process have a distinc- 
tive form (Fig. 3): 

Input: streamwise vortex, 

Output: streamwise streak. 

A streamwise vortex is an elongated region of vorticity approximately aligned 
with the x axis, and a streamwise streak is an elongated region of high or low 
velocity (relative to the mean flow) approximately in the x direction. Streamwise 
vortices and streaks are persistent features in laboratory experiments involving 
all kinds of internal and boundary layer shear flows (25). Physically, they are 
not hard to explain: in a shear flow, a small perturbation in the form of a 
streamwise vortex may move fluid from a region of higher to lower x velocity, 
or vice versa, where it will appear as a large local perturbation in the x velocity 
(26). Since these features constitute 3D perturbations of the flow field, however, 
their prevalence has been difficult to reconcile with the predictions of eigenvalue 
analysis. Non-modal analysis offers a resolution of this difficulty, for although 
streamwise streaks are not eigenmodes of the linearized flow problem, they are 
pseudomodes. 

To explain this term and our use of the words “input” and “output,” we 
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Fig. 3. Schematic illustration of a structure that appears in many shear 
flows. A low-amplitude vortex or counter-rotating pair of vortices roughly 
aligned with the flow entrains fluid from regions of high to low z velocity 
and vice versa. The entrained fluid appears as a streak of locally high or 
low streamwise velocity (not shown). This process is linear, but unrelated 
to eigenmodes. 

must define the evolution operator that is the mathematical basis of this article. 
Let u 0 = u 0 (x,y,z ) denote the vector velocity field corresponding to Poiseuille 
or Couette laminar flow (Fig. 1). Let u 0 + u(f) = u 0 (x,j/,z) + u(x,y,z,f) be the 
velocity field of a slightly perturbed flow, i.e., a nearby solution to the Navier- 
Stokes equations. (We use bold face letters to distinguish quantities that are 
functions of t .) If we take u to be infinitesimal, then it satisfies an equation 

^(f) = 

where £ is a linear operator that we call the linearized Navier-Stokes evolution 
operator. (The factor -i = -y/-L is a matter of convention.) As a measure of 
the size of solutions to this equation we define 

/ , , \i/2 

||u(t)|| = (y \u(x,y,z,t)\ 2 dxdydz) , 

which we call the energy norm since its square corresponds to energy. 

Suppose now that the linearized fluid system is driven by a signal of the 
form v(t) = e~ iu,i v for some frequency and function v = v(x,y,z): 

— (t) = -iCu{t) + e~ lult v, -oo<t<oc. 
dt 

It is easily verified that the response will be 

u(<) = u — (ljI — C) 1 t% 
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where / denotes the identity operator. Thus the operator (ujI — £) _1 , known as 
the resolvent of £, transforms “inputs” v to the linearized fluid flow at frequency 
u> into corresponding “outputs” u (10, 11). The degree of amplification that may 
occur in the process is equal to the operator norm 

I*/ t — in ll^ll 

II {uI-C) = sup—. 

«*o M 

Since an arbitrary time-dependent perturbation of the laminar flow can be re- 
duced to an integral over real frequencies by Fourier analysis, values of lv on the 
real axis are of particular interest. The maximum possible amplification over 
all real frequencies is sup^^ ||(u;/ — £)~* ||, and this is the quantity plotted in 
Fig. 2. 

An eigenvalue of £ is a number such that Cu = wu for some corre- 

sponding eigenfunction u. Equivalently, it is a number lo with the property that 
inputs with frequency uo can be amplified unboundedly: ||(u;/ — £) -1 || = oo. 
More generally, for any e > 0, an e-pseudo-eigenvalue of £ is a number u 
such that ||(u:/ — £) -1 || > e , and a corresponding e-pseudo-eigenfunction or 
e-pseudomode is any function u with ||£u— ura|| < e||u||. If e is very small, then 
an e-pseudomode u may be excited to a substantial amplitude by a very small 
input, possibly including noise in an experimental apparatus. At R = 5000, for 
example, a streamwise streak is an e-pseudomode of the linearized Poiseuille 
flow problem for e ss 1.2 x 10 -5 (Fig. 2) and thus can be excited by a streamwise 
vortex five orders of magnitude weaker in amplitude. 

The set of e-pseudo-eigenvalues of an operator is the e-pseudospectrum (27): 

A,(£) = {weC: ||(w/-£)- 1 ||> e - 1 }. 

The pseudospectra { A e (£)} form a nested family of sets in the complex plane, 
with A 0 (£) equal to the spectrum A(£). If £ is normal, A f (£) is the set of 
all points at distance < e from A(£), but in the non-normal case it may be 
much larger. The physical significance of pseudospectra can be described as 
follows. If a system is governed by a linear operator that is normal or close to 
normal — familiar examples include musical instruments, oscillating structures, 
and molecules as described by quantum mechanics — then the frequencies at 
which it resonates strongly are determined by its spectrum. The resonances of 
a system that is far from normal, however, can only be understood by examining 
its pseudospectra. 
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Spectra and pseudospectra 

Figs. 4 and 5, the centerpieces of this article, depict spectra and pseu- 
dospectra for Couette flow at R = 350, 3500 and Poiseuille flow at R = 1000, 
10000. These Reynolds numbers roughly span the range in each case from pos- 
sible turbulence (in some experiments) to unavoidable turbulence (even in ex- 
periments under the most carefully controlled conditions). These and our other 
figures are the results of numerical computations. As a practical matter we first 
Fourier transform in x and 2 , reducing the calculation to one space dimension 
(y) and two real parameters a, ft (wavenumbers in x,z). The generation of each 
figure then requires a minimization in the a-ft plane (28). 

Spectra of £ a p for particular values of a and ft have been published by var- 
ious authors (5, 29), and pseudospectra of C a p will appear in (12, 16). Neither 
spectra nor pseudospectra of C itself have been plotted previously. 

Fig. 4 reveals that in the case of Couette flow, the spectrum of £ is a contin- 
uous region contained in the lower half-plane. (For each a, ft pair the spectrum 
is discrete, but the union over a and ft is a continuum, shown in the figure as 
a shaded region.) This is true for all 7?, corresponding to the unconditional 
stability of Couette flow according to eigenvalue analysis. The spectrum comes 
closest to the upper half-plane along the imaginary axis, at a distance 0(ii -1 ) 
as R — » oo. The pseudospectra of £, by contrast, protrude significantly into the 
upper half-plane, indicating that evolution process governed by this operator 
will feature prominent non-modal effects. 

One cannot see in Fig. 4 what values of ct and ft correspond to various 
points in the spectrum and pseudospectra. The situation is roughly as follows. 
The upper boundary of the spectrum corresponds to ft — 0 and to values of a 
that grow in proportion to Rea;. In the pseudospectra, the values of a are 
comparable but the values of ft are 0(1). This difference in values of ft reflects 
the fact that although the dominant eigenmodes of this problem are 2D, the 
dominant pseudomodes are 3D. 

Fig. 5 reveals the more complicated situation for Poiseuille flow. The pseu- 
dospectra are qualitatively similar, again extending well into the upper half- 
plane (30). The boundary of the spectrum, however, is now determined by two 
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Fig. 4. Spectrum and e-pseudospectra in the complex w-plane of the lin- 
earized Navier- Stokes evolution operator for Couette flow at R — 350 and 
3500. The spectrum is the shaded region, and the solid curves, from outer to 
inner, are the boundaries of the e-pseudospectra with e= 10“ 2 ,10~ 2 * 5 , 10“ 3 , 
10" 3,5 . The spectrum lies in the lower half-plane, but the pseudospectra ex- 
tend significantly into the upper half-plane, revealing the non-normality of 
this operator. Note that the real and imaginary axes are scaled differently. 
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distinct components, each with (3 — 0. One of these looks qualitatively like the 
boundary of the spectrum for Couette flow and corresponds to perturbations 
of the flow field that are symmetric about the centerline. The other is very 
different, consisting of bumps on either side of the origin corresponding to an- 
tisymmetric perturbations of the flow. For lower values of /?, the spectrum is 
contained in the lower half-plane, but at R = 5772 the bumps cross the real 
axis. These bumps represent the mode that has received most of the attention 
in the literature, which is sometimes known as a Tollmien-Schlichting or TS 
wave. However, their negligible effect on the nearby pseudospectral contours 
reveals that this mode is oriented at a sizable angle to the remainder of the 
system (31). The pseudomodes associated with upJsO — streamwise streaks with 
a^O, (3^0 — are physically more important. 

We can summarize the implications of Figs. 4 and 5 by noting that whereas 
spectral analysis suggests that plane Couette and Poiseuille flows are fundamen- 
tally different, since one has unstable eigenmodes for certain values of R and the 
other does not, pseudospectral analysis suggests that they are fundamentally 
similar, as is observed in experiments. 


Resonance curves 

Figs. 6 and 7 represent slices of Figs. 4 and 5 along the real axis. In 
each figure ||(u>7 — C)~^ || is plotted as a function of w G K for several values 
of R, indicating the factors by which linearized Couette and Poiseuille flows are 
capable of amplifying time-dependent perturbations at various real frequencies. 
When the amplification is large we speak of the phenomenon of resonance, 
although a more careful term might be “pseudo-resonance.” 

For Couette flow, no significant resonance is present except near cj = 0. 
There, the amplification may be powerful indeed, reaching a value 0(R?) al- 
though the spectrum lies at a distance 0{R~ l ) below the real axis. For Poiseuille 
flow, the behavior for u ~ 0 is similar, but for R > 5772 an additional resonance 
for larger u appears due to the unstable eigenmode. In principle, this mode 
would dominate if one could set up an experiment involving a sufficiently long 
channel and a sufficiently smooth flow. Under conditions where this can happen, 
however, low-frequency inputs are already subject to non-modal amplification 
by a factor of order 10 5 . 
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bations to the laminar flow field at frequency u 



Fig. 7. Like Fig. 6 but for Poiseuille flow. For #>5772 infinite amplifica- 
tion is possible in principle, but of doubtful significance in practice. 
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Transient energy growth 

Another aspect of linear, non-modal instability of fluid flows is the transient 
amplitude or energy growth that may evolve from certain initial conditions. This 
is the phenomenon that has been studied in the recent papers by Gustavsson, 
Butler and Farrell, and Reddy and Henningson, where many more physical 
details can be found. 

Mathematically, consider the initial-value problem 

~^{t) = —iCu(t), t > 0, ii(0) = u. 

The solution can be written u(f) = exp(— itC)v, where exp(— itC) is the operator 
exponential. The factor by which such solutions can grow in time t is 

l|exp(— i/£)|| = supMM, 

«*<> INI 

and the energy growth is the square of this quantity. If C were a normal operator 
with spectrum in the lower half-plane, we would have ||exp(— itC)\\ < 1 for all 
t > 0. In actuality, Gustavsson and Butler and Farrell discovered that the growth 
is || exp(— it >C) || = O(R), making the energy growth 0(R 2 ). 

Fig. 8 illustrates this phenomenon for Couette flow at various Reynolds 
numbers. The agreement of the curves for finite R and R — oo can be explained 
as follows: the stream wise vortex/streak mechanism of energy growth is inviscid, 
and operates on a time scale O(R) before being shut off by the effects of viscosity. 
Such behavior is physically straightforward, appearing complicated only when 
interpreted in the basis of eigenmodes. 

Fig. 9 gives further information for Couette flow with R = 350; other 
Reynolds numbers are qualitatively similar. The restriction to x-independent 
flows (or = 0) curtails the transient growth somewhat for small t but has no ef- 
fect for larger t, when the maximal growth is achieved with a = 0 anyway (i.e., 
a purely streamwise structure). The restriction to ^-independent flows (/? = 0), 
however, has a huge effect, shutting off most of the transient growth except for 
f « 0 (32). 

The 0{R 2 ) resonances of Figs. 2,6,7 are approximately equal to the inte- 
grals under the curves in Figs. 8 and 9. Mathematically, this can be seen from 
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Fig. 8. Non-modal transient energy growth for Couette flow. For any 
finite R , the effect of viscosity shuts off the growth on a time scale 0(R). 



Fig. 9. Transient growth for Couette flow at R = 350. Little growth occurs 
when only 2D perturbations are allowed. 
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100 1000 5000 10,000 25,000 50,000 



Fig. 10. Analogue of Fig. 2, showing transient energy growth factors 
sup 1>0 || exp(— it£)|| 2 instead of resonant amplification factors sup^^ [|(u>/— 
C)~ l || . For example, transient energy growth by a factor of 1000 is possible 
for all R > 2260. If C were normal, the growth factors would be exactly 1 
outside the black region. 


the formula £ _1 = i / 0 °°exp(— it£)dt, which implies ||£ -1 || < / 0 °° ||exp(— itC)\\ dt; 
in practice the inequality is typically within a factor of 2 of equality. We can 
interpret the 0(R 2 ) result physically by noting that the resonant amplification 
is due to a combination of two effects: one (normal) factor 0(R ) representing 
the time scale over which input energy can accumulate before it eventually de- 
cays, and another (non-normal) factor O(R) representing transient amplitude 
growth. 

Analogues of Figs. 8 and 9 for Poiseuille flow look qualitatively the same; 
the differences for R > 5772 appear only for large t , off the scale of our plots. 
Fig. 10 is an analogue of Fig. 2 for the Poiseuille case that shows transient 
energy growth factors instead of resonance amplitudes (33). 

A lower bound for transient growth factors can be derived from the geom- 
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etry of the pseudospectra in the complex upper half-plane. One can show: 

sup||exp(— itC)\\ > sup e _1 a e (£), 

«>0 e >0 

where cr e (C) = sup weAt (£) Iraw is the e-pseudospeclral ordinate of C (11, 16). 
In words, if the pseudospectra protrude far into the upper half-plane, then sub- 
stantial transient energy growth must be possible. The horizontal line segment 
in Fig. 9 marks this lower bound, which falls short of the true maximum by 
only about a factor of 1.4 in amplitude or 2.0 in energy, numbers typical of both 
Couette and Poiseuille flows for a wide range of R (3, 34). 


Physically interesting pseudomodes 

Instead of operator norms, it is customary in the literature to investigate 
the evolution of particular solutions chosen for one physical reason or another. 
The amplitude history of the solution with initial flow field u(0) = v is given by 
||exp(— itC)v\\, and the upper envelope of all such curves, corresponding to all 
initial functions v , is the operator norm || exp(— if£)|| that we have discussed. 

Three choices of v with particular claims to attention can be identified 
with the aid of the singular value decomposition (SVD), the natural tool for all 
kinds of questions involving norms and extrema of nonsymmetric matrices and 
operators (35). Fig. 11 displays the curves ||exp(— ?/£)u|| corresponding to the 
principal right singular vectors of the following operators: 

iq : exp(— iO + C) (maximal initial growth rate), 

v 2 : exp (—it 0 ptC) (maximal total growth), 

v 3 : C~ l (maximal resonance). 

The function tq is the perturbation with maximal growth rate at t = 0, a function 
that has been studied by Serrin, Busse, Joseph, Lumley, and others (2, 8, 15, 
36). (The notation exp(— i0 + £) indicates the limit of the SVD of exp(— itC) 
as t — > 0. A cleaner and computationally more appropriate characterization of 
the function Uj is that it is the principal eigenfunction of (£ — C*)/2i, with 
growth rate at t = 0 equal to the corresponding eigenvalue, or equivalently, to 
the distance that the numerical range of C protrudes into the upper half-plane.) 
The function v 2 is the Butler-Farrell “optimal” that achieves maximal total 
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Fig. 11. Initial growth and eventual decay resulting from three particular 
initial perturbations v for Couette flow at R — 350. The dashed curve is 
the operator norm of Fig. 9. 

growth at some time t — t op t. The function u 3 is the one that excites a maximal 
resonant response in the sense of Figs. 2,6,7. 

Fig. 11 reveals that although the physical ideas behind the functions Uj, 
v 2 , are different, v 2 and i> 3 at least lead to comparable energy growth. These 
functions are also similar in structure: both have approximately the form of 
stream wise vortices that evolve into streamwise streaks (corresponding to left 
rather than right singular vectors of the same operators). An analogous figure 
for Poiseuille flow looks roughly similar, but the gap between v 2 and v 3 is smaller 
and the gap between these and is larger. 
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The inviscid limit 

The relationship between the stability theories for viscous and inviscid flows 
has been an awkward one (21). It is natural to expect that the mechanism of 
instability should be essentially inviscid, so that it should be possible to get at 
the essence of the matter by setting R = oo, simplifying the equations greatly. 
Yet the results of inviscid stability analysis have generally matched neither the 
results of viscous analysis nor experimental observations. A conspicuous exam- 
ple is Rayleigh’s theorem of 1880 (37), which asserts that unless the velocity 
profile contains an inflection point, an inviscid parallel shear flow is always sta- 
ble, a conclusion plainly at odds with the observed behavior of Poiseuille and 
Couette flows at large finite Reynolds numbers. 

We believe that much of the difficulty with the limit R — > oo has been due 
to misunderstanding of the significance of eigenvalues. Rayleigh’s theorem, like 
Squire’s theorem, is actually a statement about the existence of eigenvalues in 
the upper half-plane, and like Squire’s theorem, it is mathematically correct but 
misleading. If one defines stability by energy growth instead of by eigenvalues, 
one finds that every planar shear flow with R = oo is unstable, regardless of the 
velocity profile. This observation is due to Ellingsen and Palm (38), with roots 
going back to Kelvin, and related results have been established by many others 
(23, 39). The example of Couette flow with R — oo and a = 0 is particularly 
simple, for the equations reduce to the 2x2 block matrix problem 



where v and f) represent the y components of velocity and vorticity (functions 
of y, after Fourier transformation in x and z). The matrix L in this equation 
is defective (i.e., nondiagonalizable), with ||exp(— itL)\\ ~Ct as t — » oo for some 
constant C = C(&) < 1 (after transformation to the energy norm). Thus inviscid 
Couette flow is linearly unstable, even though there is no exponentially growing 
eigenmode, which explains the dashed line in Fig. 8. 

This “algebraic” instability of inviscid flows is well known among hydro- 
dynamicists. What has been missing is a fully satisfactory link to the case of 
viscous flows, where the evolution operator is typically nondefective and eigen- 
value analysis has been accepted as a test of stability. The point that has not 
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always been appreciated is that although a flow may be stable in the sense of 
bounded response to infinitesimal perturbations for each finite /?, it does not 
follow that it satisfies any bound uniformly as R — * oo. In the absence of a 
uniform bound, the notion of stability for large finite Reynolds numbers has 
little meaning (40). 

Our computations indicate that by most measures, the limit R — * oo is a 
smooth one. Table 1 summarizes the behavior of a number of quantities in this 
limit and may serve as a reference point for most of our figures (41). 



Couette 

a 

P 

Poiseuille 

a 

p 

distance of spectrum from real axis 

(72/2.47)" 1 

0 

0 

(fl/2.47)- 1 

0 

0 

max. resonance sup^ e ]^ ||(u?/ — £) _1 || 

(72/8. 12) 2 

0 

1.18 

(72/17. 4) 2 

0 

1.62 

transient growth sup i>0 1| exp(— itC)\\ 2 

(72/29. 1) 2 

35.7 /R 

1.60 

(72/71. 5) 2 

0 

2.04 

optimal time £ 0 pt 

72/ 8.52 



72/13.2 

w 

n 

lower bound based on pseudospectra 

(T2/42.6) 2 

0 

1.62 

(72/103. ) 2 

0 

2.04 

transient growth (a = 0) 

(T2/29.3) 2 

0 

1.66 

(72/71. 5) 2 

0 

2.04 


Table 1 . Leading-order behavior of various quantities as R—> oo. These 
formulas are accurate to at least 1% for all R > 100. The results for 
Poiseuille flow pertain to the highly non-normal part of the problem, ig- 
noring the R > 5772 mode (TS wave). 


Destabilizing perturbations and secondary instability 

There is an alternative, equivalent definition of the pseudospectra of an 
operator £ (27, 42): 

A e (£) = closure {u;€ C : u € A (£ + £) for some £ with ||£||<e}. 

In words, the e-pseudospectrum of C is the union of the spectra of all perturbed 
operators £ + £ with ||£|| < e, together with any limit points of this set (a tech- 
nicality of little importance). It is interesting to reconsider the implications 
of Figs. 4 and 5 in the light of this alternative characterization. The pseu- 
dospectral boundary contours in these figures imply that although the flows in 
question are eigenvalue stable (for R < 5772 in the Poiseuille case), exceedingly 
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Fig. 12. An eigenvalue stable shear flow at Reynolds number R can be 
destabilized by a linear operator perturbation Z of norm = 0(R ~ 2 ); 
the dependence is so nearly quadratic that the curves appear straight to 
plotting accuracy. The inverse e~j n = 0(R 2 ) is equal to the maximal res- 
onant amplification plotted in Fig. 2. The circles mark the approximate 
Reynolds numbers (350, 1000) at which Couette and Poiseuille flows are 
observed to undergo transition to turbulence. 

small perturbations to the evolution operator, of order 0(R~ 2 ), suffice to make 
them eigenvalue unstable. The norm of the minimal destabilizing perturba- 
tion is e min = (sup w€R ||(u;/ — £) -1 [|) — 1 , i.e., the inverse of the quantity plotted 
in Fig. 2 and labeled maximal resonance in Table 1. This observation pro- 
vides a new interpretation of Fig. 2. For example, although Poiseuille flow with 
R = 5000 is eigenvalue stable, there exists a perturbation S of norm 1.21 x 10 5 
that renders it unstable. See Fig. 12. 

This raises the question of the physical meaning of operator perturbations 
and what relevance they may have to flow instabilities observed in the labora- 
tory. (We emphasize that £ is a perturbation of £, not of the flow field u 0 .) 
The minimal destabilizing perturbation £ with ||£|| = e min is easily character- 
ized: it is a rank-1 operator of the form avu*, where a,v,u are the principal 
singular value and right and left singular vectors of £ -1 , respectively. This op- 
erator transforms streainwise streaks into streamwise vortices, thereby closing 
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the loop so that the transient growth of Figs. 8-9 can feed back upon itself to 
become modal growth for all t. Of course there is no reason to expect such 
a perfectly contrived perturbation to arise under natural conditions. What is 
suggestive is that even random perturbations of C have qualitatively the same 
effect. At R = 5000, for example, approximately 90% of all random perturba- 
tions £ to C with HFII = 10 -3 render a Poiseuille flow eigenvalue unstable, where 
by a random perturbation we mean an elementwise random perturbation of the 
20x20 matrix that represents the projection of C onto the space spanned by its 
20 dominant eigenmodes (with fixed a = 0 and j3 = 1.62). Though it is a large 
step from such matrices to shear flows, the possibility is clearly suggested that 
seemingly negligible perturbations may destabilize a flow even in the narrow 
eigenvalue sense. One can imagine that imperfections in the boundary walls of 
a laboratory apparatus, for example, might have this effect. We hope that fu- 
ture investigations will reveal whether this idea of perturbed evolution operators 
provides an explanation of instability in some flows. 

A different prospective application of the idea of destabilizing operator per- 
turbations may be to the theory of “secondary instability” as an explanation of 
subcritical transition to turbulence (8, 9). This theory is founded on the discov- 
ery that when certain laminar shear flows are perturbed by certain physically 
motivated waves of large amplitude, the resulting problem is eigenvalue unsta- 
ble, with eigenvalues high enough in the upper half-plane to match observed 
time scales for transition. The great sensitivity of the spectrum of C, to pertur- 
bations, however, raises the question of whether similar eigenvalue instabilities 
might also arise with the use of more or less arbitrary perturbations, not just 
the ones that have been considered based on physical arguments. 
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Nonlinear bootstrapping and transition to turbulence 

Kelvin wrote in 1887 (19): 

It seems probable, almost certain indeed, that... the steady motion is stable for 
any viscosity, however small; and that the practical unsteadiness pointed out by 
Stokes forty-four years ago, and so admirably investigated experimentally five 
or six years ago by Osborne Reynolds, is to be explained by limits of stability 
becoming narrower and narrower the smaller is the viscosity. 

This view of instability is still the standard one, but to this day, it has never been 
confirmed in detail. In this final section we would like to speculate about what 
the eventual confirmation may look like — in other words, about how nonlinear 
and linear mechanisms interact to bring about transition to turbulence. 

Consider the 2x2 nonlinear model problem 

d i = Au + \HBu, A=(- R 0 ' B=(" ■;). 

where R is a large parameter. The linear term, involving the non-normal ma- 
trix A , amplifies energy transiently. The nonlinear term, involving the skew- 
symmetric matrix Z?, rotates energy in the Uyii 2 plane but does not directly 
create or destroy energy, since it acts orthogonally to the motion (43). Thus 
we have linear amplification coupled with energy-neutral nonlinear mixing, a 
situation that holds also for the equations of fluid mechanics (44). 

Fig. 13 shows norms ||u(f)|| for solutions starting from eight different initial 
vectors of the form u(0) = (0,const)^, with R = 25. For ||u(0)|| < 10 , the 

curves are approximately translates of one another on this log scale, indicating 
that the evolution is effectively linear. At ||u(0)|| =4 x 10“ 4 , the nonlinearity 
has begun to have a pronounced effect. At ||u(0)|| = 5 x 10~ 4 the threshold has 
been crossed, and for this and higher amplitudes, the curves do not decay to 
zero but blow up to a critical point of amplitude ~ 1. 

Fig. 13 reveals a remarkable phenomenon: the amplitude growth is far 
greater than that of the purely linear problem dujdt = Au. Table 2 lists the 
linear growth factor M ~ Rj\ and the threshold amplitude e (for this particular 
choice of initial vector) for four values of R. Evidently e. is of order M 3 , not 
M -1 . This “bootstrapping” effect can be explained as follows. Suppose the 
solution at t = 0 consists of a vector of amplitude e in a direction that excites 
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Fig. 13. ||u(f)|| vs. t for solutions to a nonlinear 2x2 model problem with 
various initial amplitudes. The threshold amplitude is ||u(0)|| « 4.22 X 10 -5 . 


R 

linear 

growth 

factor 

M 

threshold 

amplitude 

€ 

12.5 

3.18 

3.41 x 10~ 3 

25 

6.28 

4.22 xlO" 4 

50 

12.5 

5.27 xlO" 5 

100 

25.0 

6.58 x 10~ 6 


ratio 


8.1 


8.0 


8.0 


Table 2. The “bootstrapping” of the linear growth factor M = 0(R ) 
causes the threshold amplitude to be 0(R~ 3 ). 


growth of the linear problem du/dt = Au (the principal right singular vector of 
>4). At a later time, of order i?, the solution has grown to order Re by the linear 
growth mechanism but moved into a direction that no longer excites growth (the 
corresponding left singular vector). Meanwhile, however, the nonlinear term has 
had the effect of transferring some of this energy back to the original direction, 
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with amplitude R(Re) 2 = R 3 t 2 since the nonlinearity is quadratic and the time 
scale is O(R). If R 3 c 2 is of order less than e, the process is not self-sustaining 
and the energy dies away. On the other hand if R 3 e 2 is of order greater than e, 
then there is more energy than was present at the start and feedback occurs, 
leading to self-sustained amplitude growth. Thus the threshold amplitude is 

e=0(R~ 3 ). 

A similar experiment shows that if the same nonlinear equation is driven 
by a forcing oscillation instead of an initial vector u(0), the threshold 

amplitude becomes e = 0(R~ 4 ). 

It may appear at first sight that these results indicate the great power and 
importance of nonlinear effects. Yet in two senses, these energy growth scenar- 
ios are essentially linear. First, as mentioned above, the nonlinear term does 
not add energy locally but merely redistributes it. Second, the appearance of 
the bootstrapping phenomenon does not depend on the precise nature of the 
nonlinearity. Any quadratic nonlinear term that transfers energy from decaying 
to growing solution components has the potential of inducing a threshold ampli- 
tude t = 0(R~ 3 ) or 0(R ~ 4 ) with respect to initial or forcing data, respectively; 
a random perturbation, for example, will often suffice. Higher-order nonlinear- 
ities also lead to similar effects, though the exponents may be lowered, e.g. to 
e= 0(R~ 2 ) and 0(R~ 3 ) for a cubic nonlinearity. 

The Navier- Stokes equations are certainly more complicated than our 2x2 
model. One difference is that instead of 2-vectors, they act on functions with 
infinitely many degrees of freedom, most of which do not experience non-normal 
linear growth. There will always be some energy in the growing pseudomodes, 
however, and in a pipe or channel of substantial extent, random fluctuations 
can be expected to raise the energy levels in such components locally well about 
the statistical average (45, 46). Another difference is that the nonlinear in- 
teractions in the Navier-Stokes equations act across different wave numbers a 
and j3 (13, 14), making the exponent of 2 in our model perhaps unrealistically 
low. Despite these qualifications, we conjecture that the transition to turbu- 
lence of eigenvalue-stable shear flows proceeds analogously to our model in that 
the destabilizing mechanism is essentially linear in the senses described above, 
with an amplitude threshold for transition of order 0(R a ) for some a < — 1. 
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Conclusion 

In this article we have described three essentially linear approaches to the 
phenomenon of instability of shear flows: (i) (pseudo-) resonance, (ii) transient 
energy growth, and (iii) destabilizing operator perturbations. These ideas are by 
no means independent. Mathematically, they are all related to the pseudospec- 
tra of the operator £, and physically, they all depend on the same mechanisms 
of extraction of energy from the mean flow by streamwise vortices and other 
structures. One should not expect that any one of these ideas will necessarily 
prove to be the “right” one, but rather, that each of them may prove appropriate 
to a particular class of experiments. For example, it is reasonable to speculate 
that in different situations transition to turbulence may be primarily triggered 
in physically quite different ways by (i) laboratory vibrations, (ii) initial or in- 
let disturbances, or ( Hi ) deviations of the pipe or channel geometry from the 
Poiseuille or Couette ideal. 

It is likely that in the upcoming decade, great progress will be made in elu- 
cidating these details. Numerical simulations of the 3D Navier-Stokes equations 
are becoming routine, and in the next few years such simulations will probably 
settle at last the question mentioned in the last section, for example, of how the 
threshold amplitude for transition depends on the Reynolds number in various 
flow situations. 

Besides hydrodynamic stability, there are other fields also where highly 
non-orthogonal eigenfunctions arise and eigenvalues may be misleading. Exam- 
ples in fluid mechanics include the instability of magnetic plasmas (48) and the 
formation of cyclones (49). Two more examples come from the field of numeri- 
cal analysis: the numerical instability of discretizations of differential equations 
(40) and the convergence of iterative algorithms for nonsymmetric systems of 
linear equations (50). The recurring theme in these and other applications is 
that although the long-time behavior of an evolving system may be controlled 
by nonlinearities, some important phenomena are of a short-time nature and 
essentially linear (51, 52). If the linearized problem is far from normal, eigenval- 
ues may be precisely the wrong tool for analyzing it, for eigenvalues determine 
the long-time behavior of a non-normal linear process, not the transient. We 
can expect to see further examples of alternatives to eigenvalues in the years 
ahead. 
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